ON THE ASYMPTOTIC NUMBER OF EDGE STATES FOR 
MAGNETIC SCHRODINGER OPERATORS 
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■ Abstract. We consider a Schrodinger operator (hD — A) 2 with a pos- 
itive magnetic field B = curl A in a domain Q, C R 2 . The imposing of 
Neumann boundary conditions leads to spectrum below h'mi B. This is 
a boundary effect and it is related to the existence of edge states of the 
system. 

£ — 1 We show that the number of these eigenvalues, in the semi-classical 

limit h — > 0, is governed by a Weyl-type law and that it involves a 
symbol on dQ. In the particular case of a constant magnetic field, the 
curvature plays a major role. 

O 

■ 1. Introduction and main results 
O , 

\& , 1.1. Introduction. In this paper we consider a magnetic Schrodinger op- 

erator 

(1.1) P h := (KD - A) 2 in L 2 (n) 

where D = — iV and £1 C M? has a smooth and compact boundary dfl 
on which we impose Neumann boundary conditions. The 'magnetic field' 
curl A = B is assumed to be smooth, positive and 'not too large on the 
boundary' (see Ql.tjj) below for the precise assumption). As we will see, the 
choice of Neumann boundary conditions implies the existence of eigenvalues 
below h'mfB. Our goal is to determine for bo < inf B the asymptotic 
H ' behavior of the number 



(1.2) N(hb ,P h ) 

of eigenvalues of below hbo in the semi-classical limit h — ► 0. The semi- 
classical limit is of course equivalent to the limit of a strong magnetic field. 

The operator (|1.1|) has recently received a lot of attention in connection 
with the Ginzburg-Landau theory of superconductivity. Indeed, the lowest 
eigenvalue of (|l.lj) plays an important role in the description of the boundary 
nucleation of superconductity close to the third critical field, see fLuPa , 
|HeMoj . [FoHc2 and references therein. The methods developed in these 
papers will also be important for our analysis. They show in particular 
that the boundary has an effect similar to that of a potential well and that 
eigenfunctions of Ph corresponding to eigenvalues below h inf B are strongly 
localized near the boundary. Hence they correspond to 'edge states'. 
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We emphasize that the 'energy level' in (|1.2j) depends on h. If instead the 
energy A > is fixed, it is well-known that 

N(\,P h )~h- 2 \^ 

for bounded f2. In particular, the leading term is independent of the mag- 
netic field. Recalling that the distance between two Landau levels is pro- 
portional to h it is natural to expect effects of the magnetic field to appear 
when A = boh. We note also that (|1.2j) can be rewritten as 

N(hb , P h ) = N(b , {h x ' 2l D - h- l ' 2 Af). 

This suggests that the effective semi- classical parameter of the problem is 
/i 1 / 2 , which tends to zero simultaneously as the effective magnetic field 
h~ l l 2 B tends to infinity. 

The first results on spectral asymptotics of the operator (jl.lj) were ob- 
tained by Colin de Verdiere |CVj and Tamura [Taj . The former author 
considers a slightly different problem, namely the behavior of N(X,Pi) as 
A — ► oo. Moreover, Pi is defined there as the Dirichlet realization of (jLlj) . 
The methods of |CV| (see also (Tr]) allow to prove that under suitable as- 
sumptions on A one has for all bo > 



oo 

(1.3) N^P^-h- 1 —^ B(x)dx. 

Z7T ^[Jixan: (2n-l)B(aO<M 

This formula is valid both for the Dirichlet and for the Neumann realization. 
However, while the Dirichlet realization has no spectrum below Ziinf B, this 
is no longer true for the Neumann realization. In this case (|1.3j) provides 
the estimate 

N{hb ,P h ) = o(h- 1 ), b <MB. 

Our goal is to improve upon this and to determine the precise behavior of 
N(hbo,Ph) as h — > 0. So loosely speaking we are interested in the leading 
term of the spectral asymptotics whose first term vanishes. 

For further results on spectral asymptotics of magnetic Schrodinger op- 
erators we refer to LiSoYnj, [E], |DiE,aj and references therein. 

1.2. Main result. Let us state the precise assumptions on Q and A under 
which we shall work henceforth. We assume that Q C M 2 is an open domain 
with dQ, £ C 4 . Moreover, for the sake of simplicity, we assume that dQ, is 
bounded and connected. Note that we do not assume that is bounded or 
simply connected. We consider a 'magnetic vector potential' A £ C 2 (S7,M 2 ) 
and introduce the 'magnetic field' B := curl A and the quantities 

b := inf B(x), b' := inf B{x). 

The operator ()1,1|) is defined via the quadratic form 

q h [u) ■= / \(KD - A)u\ 2 dx 
Jn 
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with domain 

V[q h ] := {u e L 2 (n) : (WD - A)u G L 2 (fi, C 2 )}. 

If the essential spectrum of Ph is non-empty, which may happen if Vt is 
unbounded, we use the well-known inequality 

(1.4) hb f \u\ 2 dx< [ | (WD- A)u\ 2 dx, u G Cq°(Q), 

in in 

and a 'magnetic version' of Persson's lemma to conclude that inf a ess (Ph) > 
hb. Hence in any case for A < hb the spectrum of Ph in the interval [0, A) 
consists of finitely many eigenvalues of finite multiplicities, and we denote 
their total number (taking multiplicities into account) by N(X,Ph). 

In order to state our main results we need some notation. For (£lwe 
denote by ^(£) the lowest eigenvalue of the operator 

-^ + (£ + i) 2 inL 2 (M+) 

with a Neumann boundary condition at the origin. Then (see Subsection l2.ll 
for more details) the minimum 

6 := inf//(0 



is attained at a unique £o an d one has £o £ ( — 1)0) and 

(1.5) Cj := //'(£o)/6M > 0. 
Throughout the following we shall assume that 

(1.6) < Q b' < b. 

Since < Oo < 1 (numerically @q = 0.59 . . . ) this is in particular true in 
the important special case of a constant magnetic field. 
Our first main result is the following. 

Theorem 1.1. Let Q b' <b <b. Then 

(1.7) Ymxh 1 ' 2 N{hbQ,P h ) = ^- ! f B(x) 1/2 ds(x)d£. 

2?r J J{(x,OednxR: B(x)^)<b } 

We emphasize that (|1,7|) has a Weyl-type form, involving the symbol 
B(x)[i(!;) on the co-tangent bundle of the one-dimensional manifold dQ. 
The essentially one-dimensional nature of the asymptotics is also reflected 
in the fact that the effective semi-classical parameter h 1 ! 2 appears with 
the power —1 in the asymptotics. This should be compared with for 
6o > b. There both the integral term and the power of h 1 ! 2 reflect the 
iwo-dimensional nature of the bulk states. 

Note that TheoremHTIimphes that N(hQ b', Ph) = o(/i -1 / 2 ). It is natural 
to ask whether it is possible to obtain the correct asymptotics. This will 
probably involve the geometry of the set {x £ dtt : B{x) = b'}. Here we 
give an answer in the particular case where the magnetic field B is constant. 
Indeed, we do not only give the asymptotics of N(hbo,Ph) for bo = QqB 
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but we allow fro to vary with h on the scale h 1 / 2 . The result will involve the 
curvature k : R/|<9f2|Z — ► R, see Subsection 12.21 for our notation. We follow 
the usual convention that k > if O is convex. We will prove 

Theorem 1.2. Assume that B is constant and let kq E R. Then 
(1.8) 

nl/4 r\d€l\ 

lim h 1/A N(hQ B + h 3/2 dB 1/2 K , P h ) = == / (k{s) + Hi )+ ds. 

There is an important difference between Theorems 11.11 and 11.21 If O is a 
bounded domain we can consider both the interior problem $7 = O and the 
exterior problem O = R 2 \ O. Now if b$ > @Qb' it follows from (|1.7j) that 
the leading order terms of N(hbo, P^) for both problems coincide (provided 
the magnetic fields coincide on the boundary). This is no longer true if B 
is constant and 6o = ©o-B- Indeed, the asymptotics are 'complementary' in 
the following sense: for the interior problem N(hOoB, P^) is, up to leading 
order, determined by the convex part of the boundary (where the curvature 
of O is positive) and for the exterior problem N(h@oB, P^) is determined by 
the concave part (where the curvature of the obstacle O is negative). This 
observation is in the same spirit as the considerations on spectral duality in 
[HoSmj . 

Remark 1.3. Let us mention an immediate generalization of Theorems ll.il 
11.21 and their proofs. If dO, has finitely many connected components and 
if one imposes on each of them either Dirichlet or Neumann conditions, 
then a formula similar to (|1.7|) holds, but the integration is restricted to 
the Neumann components. Moreover, it is enough that only the boundary 
of the Neumann components is C 4 , and also that the assumption A G C 2 
holds only in a neighborhood of those components. 

Moreover, it would be desirable to remove the assumption of smoothness 
of the boundary. If <90 is piecewise smooth one can probably use the methods 
from |Bonj . 

1.3. Outline of the paper. The proofs of Theorems 1 1 . 1 1 and 1 1 . 21 are tech- 
nical but the main idea is rather simple. To show (|1.7|) we localize the 
problem, following HeMo , to a tubular neighborhood of the boundary of 
normal size h 3 ^ 8 and cut this into boxes of tangential size /i 3 / 8 . In each of 
these boxes we approximate the magnetic field by a constant one, see Sub- 
section 13.31 This reduces the problem to the analysis of the model problem 
of an operator with constant field in a rectangle with Neumann boundary 
conditions on one edge and Dirichlet boundary conditions elsewhere. For 
this operator we cannot separate variables but it turns out (Subsection 13. 2j) 
that its spectral counting function is fairly close to the one of the operator 
on an infinite half-cylinder. The latter problem is treated in Subsection 13. II 
The reduction to the model problem and a careful estimate of the remainder 
is achieved in Subsection 13.41 In Subsection 13.51 we complete the proof of 
Theorem ll.il We will even obtain a remainder estimate. 
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The proof of Theorem 11.21 follows a similar pattern and we will be rather 
succinct there. The analysis of the model problem is, however, significantly 
more difficult, see Subsection 14.11 

On a technical level we mention that classical Dirichlet-Neumann brack- 
eting is not possible in our situation, since additional Neumann boundary 
conditions would produce too many additional eigenvalues. This difficulty 
is overcome in |CV] by a localization technique based on the IMS formula. 
We emphasize once more that the paper |CVj concerns Dirichlet boundary 
conditions, so that the boundary effects of our Theorems 11.11 and 11.21 were 
not present there. 

1.4. Acknowledgements. The author wishes to thank Prof. B. Helffer 
for the invitation to Orsay and numerous fruitful discussions. He is also 
grateful to S. Fournais and A. Hansson for useful remarks. Financial sup- 
port through the ESF Scientific Programme in Spectral Theory and Partial 
Differential Equations (SPECT) as well as through the European Research 
Network "Postdoctoral Training Program in Mathematical Analysis of Large 
Quantum Systems" (Contract Number HPRN-CT-2002-00277) is gratefully 
acknowledged. 



2. Auxiliary material 

2.1. A family of ordinary differential operators. For any ( £ K we 

consider the operator 

(2-1) _^ + (e + t )2 inL 2 (R + ) 

with a Neumann boundary condition at the origin and denote its lowest 
eigenvalue by /«(£;). (Here we use the convention of FoHcl ; we note that in 
HcMo; £ is replaced by — £.) 

The dependence of /z(£) on £ has been studied in [DaHej (see also jBolj ) 
where the following facts were established: The function \i is smooth and 
satisfies lim^_ v _ co //(£) = 1, lim^ +oc ^(£) = +co. There exists a £o £ 
(—1,0) such that \x is strictly decreasing in (— oo,£ ) and strictly increasing 
in (£o, oo). Moreover, at £o it has a non-degenerate minimum and one has 

e := = £o- 

In view of these facts it is possible to introduce two inverse functions 
v+ : [0 O , oo) -> [fo,°o), v- : [©o,l) -> (-oo,£ ] satisfying 

fiov± = id, u + o fi\ Ko 0o) = id, v- o fj,\ (.-.oo^] = id. 

At some points below it will be technically convenient to extend the functions 
v ± by £o to the interval [0, Go). 

Recalling that the minimum of /i is non-degenerate one easily establishes 
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Lemma 2.1. For any < e < 1 there exists a constant C > such that for 
all <[3<f3'< 1 - e one /ios 



As a consequence one obtains the useful estimate 



(2.2) !/ + (^) - v-tf) < - v-{fi) + ICJW^P- 

Finally, we denote by the second eigenvalue of the operator (|2.1j) and 

put 

0! :=inf Ml (0. 



Numerically, one finds Bi = 2.63 ... . 1 Below we shall only use the following 
bound on 0i, the proof of which is due to B. Helffer. 

Lemma 2.2. One has the inequality 01 > 1. 

Proof. Denote by A(£) the first eigenvalue of the operator (|2.1[) with a 
Dirichlet boundary condition at the origin. Then by Sturm-Liouville theory 
> ^(0 f° r an £ £ R- Moreover, the variational principle implies that 

(2.3) X (0>^o-(-^ + (ti + t) 2 

where the operator on the RHS is defined in I^QR)- Using translation in- 
variance and the well-known result for the harmonic oscillator one finds that 
the RHS of JEU equals 1. □ 

2.2. Boundary coordinates. Here we would like to recall the definition 
of coordinates near the boundary of f2. Recall that we assume that d£l is 
connected, C 4 -smooth and of length 

C := \dVL\. 

Let 7 : R/£Z — ► dtt be a parametrization of the boundary with |7 ; | = 1 and 
let v(s) be the interior unit normal at the point j(s). The parametrization 
can be chosen such that det(7 ; , u) = 1, so that the curvature k is given by 

It is well-known that for sufficiently small to > the map : R/£Z x 

(o,to)->n, 

S(a,t) := 7 (s)+ti/(a), 
defines a diffeomorphism between M/£Z x (0, to) an d its image 

$(R/£Z x (0,i )) ={i£fi: dist(x,0) < t } =: U to . 

We now fix a constant k € K. (For the proof of Theorem 11.11 it will suffice 
to take k = 0.) Denoting 

(2.4) a(s,t) := 1 - t/e(s), (s, t) € K/£Z x (0, t ), 



^The author would like to thank V. Bonnaillie-Noel for this calculation. 
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and 

(2.5) a R (t) := 1 - tk, t G R + , 
we define for any u G L2(O <0 ) 

(2.6) u(s,i) := (a(s,t)/a«(t)) 1/2 u($(s,t)), (s,t) G x (0,i )- 

This induces a unitary operator from L2(^ ) to ^((0, £) x (0, to), a^dsdt). 
As in the case k = considered in |HeMoj (Appendix B) one finds that if 
u G £>[(j7i] and suppu C O, to then 

(2.7) 

?aM= I [\a- 2 \(hD s -A 1 )v\ 2 + \(hD t -A 2 )v\ 2 + h 2 W k \v\ 2 )a k dsdt 
Jo Jo 

[ C (K-k)\v(.,0)\ 2 ds, 



2 Jo 
where 

(k(s)-k)(k(s) + k(1-2^(s))) tK"(s) 5t 2 n'(s) 2 

Wk{ - S ^> : ~ 4a s (t) 2 a(s,t) 2 2a(s,t) s ~ 4a(s,t) 4 ' 

We do not give the expression for A = (A±, A2) T but note only that 

B(s,t) := d s A 2 (s,t) - d t A 1 (s,t) = a(s,t)B($(s,t)). 

3. Proof of Theorem 11.11 

3.1. The model operator on a half-cylinder. We fix B,S > and 

consider the operator P^' B = {h~D — BAq) 2 in L2((0, S) x M + ) with periodic 
boundary conditions at s G {0, S} and Neumann boundary conditions at 
t = 0. Here and in the sequel 

(3.1) Ao(a,t) := (-t,0) T . 

Remark 3.1. Since we impose periodic boundary conditions, we actually 
work on the half-cylinder R/S'Z x M + . Since this is a not a simply connected 
manifold, the magnetic field alone does not determine the operator (up to 
unitary equivalence), but one also needs to specify the circulation of A 
around the boundary M/SZ x {0}. 

We shall use the following notation for a self-adjoint and lower semi- 
bounded operator T. If Et(A), A C M, is the spectral measure associated 
with T, we put 

N(X,T) := dim ran £ T ((-co, A)) , A G M. 



If the spectrum of T below A is discrete, then N(X,T) coincides with the 
number of eigenvalues (counting multiplicities) below A. 
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Lemma 3.2. Let A < hB. Then 

N(\,P*> B ) = # |zn^ [v- (h-^x) ,* + (h-'B-ix))) . 

Recall that the functions v± are extended by £o to [0, ©o)- In the state- 
ment of the lemma we use the notational convention that (y- {(3) , v + (/?)) = 
(£o,£o)=0if/?<©o. 

~ S B 

Proof. By separation of variables the operator P.' in L2((0,S) x M+) is 
unitarily equivalent to the direct sum 

© (~ h2 §2 + (2vrn/i,S- 1 + Btf^j in ^ 0L 2 (R + ) 



(with Neumann boundary conditions at the origin). Applying the dilation 
r = hr x l 2 B x l 2 t we obtain the unitary equivalence 

-tfljL + ^nhS- 1 + Bt) 2 * hB (~ + (lirnh^B-^S- 1 + rf 

By the facts mentioned in Subsection 12.11 (particularly Lemma 12. 2 J) we con- 
clude that for A < /i-B 

N(X,P^ B ) = |{n 6 Z : n(2Trnh 1/2 B^ 2 S' 1 ) < /i _1 S -1 A} 

= tJ (Z n (2vr)- 1 /i- 1 / 2 5 1 / 2 5 (^(/T^^A), iz+^fl^A))) , 
which is what we claimed. □ 

"SB 

We note that the proof shows that N{hB , P h ' ) = +oo. Moreover, we 
easily deduce from Lemma 13.21 that for all bo < B 



(3.2) 



h l / 2 N(hb Q ,P^ B ) - (v + (B-%) - u-(B-%)) 



< h 1 / 2 . 



The relation Q3.2|) is in formal accordance with Theorem II . II since 



2vr v ^ v ' v " 2tt 



{(s,0: B M (0<6o} 



3.2. The model operator on a Dirichlet strip. We fix S, T, B > and 

consider the operator P^' T ' B = (KD - BA ) 2 in L 2 ((0,5) x (0,T)) with 
Neumann boundary conditions on t = and Dirichlet boundary conditions 
on the remaining part of the boundary. Recall that Aq was defined in H3.1|) . 

STB 

Our goal is to compare the eigenvalue counting function for P h ' ' with that 
ofP*> B . 

Proposition 3.3. For all B,T, S > and all A < hB one has 

N(X,P^ B )<N(X,P^ B ). 
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Moreover, there exists a C > such that for all B,T, S > 0, all 5 G (0, 5/2] 
and all A < hB 

N(X, P*' T ' B ) > ±N(\ - Ch 2 (5~ 2 + T~ 2 ),P 2 h (S - 8) ' B ). 

STB 

Proof. The extension by zero of a function in the form domain of P. ' ' lies 

~ S B 

in the form domain of P, ' , and the values of both forms coincide for such 
a function. Hence the first assertion follows immediately by the variational 
principle. 

To prove the lower bound we follow the ideas of |CVj . For any < 5 < 5/2 
we choose a smooth partition of unity on H./2(5 — <5)Z, 

(<^) 2 + (<^) 2 = 1 on R/2(5 - S)Z, 

such that 

2 

supple [0,5], supple [5 -6, 25 - S], £ |(^)f < Cl <T 2 . 

i=i 

The constant c\ > can be chosen independently of 5, 5. Similarly, for each 
T > let Vo, V>i fulfill 

(3.3) (V#) 2 + (Vf) 2 = 1 on R + 

and 
(3.4) 

l 

supp^ C [T/2,+oo), suppVf C [0,T], ^|W>f)'| 2 < c 2 T- 2 

with C2 > independent of T. Finally put Xi' T ( s ,t) '■= fii^^Tft): * = lj 2, 
and X ' T ( S ^) : = ^oX*)- 

Let u be in the form domain of Pf^ S . Then from the IMS formula 
with I s := (0, 2(5 - 5)) we get 

/ | (KD-BA )u\ 2 dx 

2 . 2 

= £ / |(^D-SAo)xJ'V^-^ElH V ^ 

,=0^ xR + 1=0 
2 , 

> V / \(KD - BA ) x fu\ 2 dx - c 3 h 2 (5- 2 + T 
i=0 Jl s xK + 



(3.5) 



<5,Ti i|2 



ST STB 

The function Xi' u belongs to the form domain of P h ' ' and, since u is 

(5 T 

periodic, the function xi u belongs to the form domain of the operator 

STB 

TssPfr ' Tg_ s in L 2 ((5— 6, 25— 5) x (0, T)), where rs_,5 denotes translation 

STB 

by 5— 5 with respect to the variable s. Of course, TssP h ' ' Tg_ 5 is unitarily 

STB 8 T 

equivalent to P h ' ' . Finally, Xo u belongs to the form domain of the 
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operator pjf 5),B := (KD - BA ) 2 in L 2 {I S x (T/2,oo)) with Dirichlet 
boundary conditions at t = T/2 and periodic boundary conditions at s G 
{0, 2(5 — 5)}. (We do not reflect the dependence on T in our notation, for 
operators with different T are indeed unitarily equivalent.) Hence we deduce 
from (|3.5[) by the variational principle 

N(X - c 3 h 2 {6~ 2 + T~ 2 ), P? S ~ 5) ' B ) < 2JV(A, P^ T ' B ) + iV(A, pgf )■ 

From an inequality similar to (|1.4|) with f2 replaced by R/2(5 — <5)Z x R + 
we infer that Pl^ 5),B > hB. Hence N(\,P£f) = if A < /iJ3 and the 
proposition is proved. □ 

3.3. Estimates near the boundary. Now we consider a general domain 
SI C M 2 such that dQ is bounded, connected and C -smooth. We will 
approximate the quadratic form locally near the boundary by a qua- 
dratic form corresponding to a constant magnetic field. For this we use the 
boundary coordinates (as well as the notation) introduced in Subsection 12. 21 
Throughout this section we will assume that k = 0. 

Let T G (0, to] and S G (0, C) (below T, S will depend on h and tend to 
as h — > 0). We are interested in u G T>[q^\ such that the corresponding v, 
defined in (|2.6jl . satisfies 

(3.6) suppv C [0,5] x [0,T]. 

First we use a gauge transformation to make the field on [0, S] x [0, T] 'almost' 
constant. Indeed, recalling that A G C 2 (ft,M?), one readily obtains 

Lemma 3.4. There exists a constant C > such that for all S G (0, £), 
S G [0,5] there exists a function (ft G C 2 ([0, 5] x [0,to]) such that 

A(s, t) - V<ft(s, t) = (-Bt + /3(s, t), 0) T , (s, t) G [0, 5] x [0, t ], 

where B := B(S,0) and for any < T < to 

(3.7) sup \/3{s,t)\ < C(S 2 + T 2 ). 
( s ,t)e[o,S]x[o,T] 

If u is supported in a small subset near the boundary the previous lemma 
allows us to express Q/Ju], up to a small error, via a quadratic form corre- 
sponding to a constant magnetic field. 

Lemma 3.5. There exists a constant C > such that for all S G (0, £), 
S G [0, S], T G (0, to], e G (0, 1] satisfying e > CT and for all u G T>\q^\ such 
that the corresponding v satisfies (|3.6|) one has 

q h [u}-\\{hD-BA )e-^ h v\\ 2 

< e\\{hB - BA )e- i ' t ' /h v\\ 2 + Ce~ l ((5 2 + T 2 ) 2 + h 2 ) || e -^[| 2 . 

Here B := B(S, 0) and (ft is the function from Lemma \S.J\ 
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Proof. We write In view of (|2.7|) we decompose 

q h [u] - HOD - BA )w\\ 2 = h + I 2 + I 3 , 

where we define 

Ji:= J C J° (a- 2 \(hD s -A 1 )v\ 2 + \(hD t -A 2 )v\ 2 ^ dsdt 
- \\(h-D-BA )w\\ 2 , 



I 2 -=h 2 [ [ ° W \v\ 2 dsdt, 
Jo Jo 



h 2 r c 

/ 3 :=-y J k\v(.,0)\ 2 ds. 

We begin with the (easier) terms I2 and I3. Since Wo is bounded, we obtain 
that for some constant c\ > 

(3.8) \I 2 \ < ci/i 2 ||w|| 2 . 

Moreover, for any s£ [0, L] one has 

r T 

\w(s,0)\ 2 = -2Re / —w(s,t)w{s,t)dt 
Jo ot 

< [ (e^Dtwis,^ 2 + e^ 1 \w(s,t)\ 2 ) dt. 





Since k is bounded one easily concludes that there is a constant C2 > such 
that 

(3.9) |I 3 | < e||/iA^|| 2 + c 2 e- 1 /i 2 |H| 2 

for any e > 0. Now we turn to the term I\. First we note that by Lemma l3.4l 
one has 

I x = \\a-\hD s + Bt- (3)w\\ 2 - \\(hD s + Bt)w\\ 2 . 
We use that for some constant c 3 > 

\a~ 2 -l\<c 3 T on [0,5*] x [0,T], 
and hence for all €2 > 

|/i| < (1 + c 3 T)\\(hD s +Bt- f3)wf - \\(hD s + Bt)w\\ 2 

< ((1 + c 3 T)(l + e 2 ) - l)\\{hD s + BtH| 2 + (1 + c 3 T)(l + e^ l )\\(iw\\ 2 . 

In particular, if e 2 > T then (1 + c 3 T)(l + 62) — 1 < 0462- Recalling ()3.7|) we 
obtain 

(3.10) < c 4 e 2 \\(hD s + Bt)w\\ 2 + c 5 (S 2 + T 2 ) 2 (l + 1 )||^|| 2 . 

The assertion now follows easily by setting 0462 = e and summing ()3.8|) . 
(1331) . (I3~TU1) . □ 
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3.4. Bracketing. Now we combine Proposition l3.3l and Lemma f3.5l in order 
to obtain a two-sided estimate of N(hbo, Ph) in terms of the spectral counting 
functions of operators on a half-cylinder with constant magnetic fields. 
For N £ N put 

(3.11) 5 := — and s n := nS, n = 0, . . . , N. 

We choose 

N= [h- 3 /% 

and note that 5 = 0(h 3 / 8 ). Then the core of Theorem 11.11 is contained in 

Proposition 3.6. Let < b. Under the above hypotheses there exists a 
constant C > such that for any < h < C~ l , 5 £ (0,5/2] and S n £ 
[s n _i, s n ], n = 1, . . . , N , one has 

1 N 

-J2 N (hbo ~ Ch 2 5- 2 ,P 2 h (S ' 5) ' B ") < N(hb ,P h ) 



2 n=l 



N 

\S+25,B n 



< N(hb + Ch 2 5~ 2 , ps+ 2S ' B ") 

n=l 

where B n := B(S n ,0). 

Proof. We begin with the proof of the lower bound. We apply Lemma 13.51 
with 

T = h 3 / 8 , e = h 1 ^ 

not only on [0,5], but on any [s n _i,s n ]. (It is evident from the proof of 
that lemma that the constants there can be chosen independently of n.) 
Moreover, recall that 5 = 0(h 3 / 8 ). It follows that, for some constant c\ > 
and for all u £ T>[q^\ which vanish on 

JV 

[i£(l: t{x) >T}u[J{x £tt: < t(x) < T, s{x) = s n }, 

n=l 

one has the estimate 

N rsu rT 

q h [u}<(l + h l ^)Y, / \{hT> - BnA^e-^vf dsdt 

+ c 1 h 5 / 4 \\e- i ^ h v\\ 2 , 
where v is defined by Q2.6|) . Hence by the variational principle 



n=l V 



We apply Proposition l3~3l with 5 £ (0,5/2]. since 

hb - erf* ^ a + r _ 2) > hbQ _ rfr2 



1 + /lV4 
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we obtain 

1 N 

N(hb ,P h ) > - Y,N(hb - c 2 h 2 5- 2 ,P 2 h (S ~ 5) > B " 



2 

n=l 

This is the desired lower bound. 

We turn to the proof of the upper bound. We choose for every 5 £ (0, 5/2] 
a partition of unity on M./CX, 

N 

Yj&i) 2 = 1 ° n R / CZ With SU PP C [s n -l ~S,S n + 5] 
i=l 

and such that 

N 

£i(^)'i a <<*r 2 . 

The constant C3 > can be chosen independent of 5, 5, iV. Moreover, let 
^(f, ipf be as in (pOJ) . (jH3J, and for T G (0,i ) put 

X £ T (x) := ^( S (x))^f(t(x)), n = 1, . . . , iV, x J T (x) := ^(t(x)). 

By means of the IMS formula we find for all u £ T>[q^\ 

N N 

n=0 n=0 

2V 



>Y,^[x & n T u]- c 4 h 2 (6- 2 + T- 2 )\\u\ 



n=0 



As in the proof of the upper bound choose T, e as in (|3.4j) . Then we obtain 
from Lemma 13.51 that for n = 1, . . . , N 

q h [x S n T u\ > (l - h 1 ' 4 ) f n [ T \(KD - B n A )e-^l h X s > T v\ 2 dsdt 

V 7 Js n -l Jo 

where u is related to u by (|2.6f) and where we write Xn T on the RHS instead 



of Xn *£• It follows that 



9/! 



N > (l - /i 1/4 ) V / \(hD- B n A )e-^ h x 5 fv\ 2 dsdt 

V 7 1 is„_, JO 



-1 J Sn-l JO 

1 r S , T 1 i,2r- 2ii i|2 

+ <MX U J ~ C 6^ d \\ u \\ ) 

and hence by the variational principle 

^ ' > ' 2r-2 



AT(A, P ft ) < £ N ( X + **£\ P? + * TA ) + iV(A + c,hH- 2 , P£ 

n=l ^ 7 



14 RUPERT L. FRANK 

where P^ Q := (KD — A) 2 in L2(Q\Qt/2) with Dirichlet boundary conditions. 
Because of the inequality (|1,4|) we have N(X, P^q) = for A < hb. Moreover, 
since 

hb + c 6 h 2 5- 2 , , ,2,-2 
_ 1/4 — < hb + c 7 h 6 , 

the desired upper bound follows from Proposition 13.31 □ 

3.5. Proof of Theorem 11.11 To complete the proof of our main result 
we now combine Proposition 13.61 with the explicit calculation in Lemma 13.21 
for the half-cylinder. Indeed, instead of (|1.7[) we prove the stronger (albeit 
probably not sharp) estimate 

(3.12) h^Nlhb^Ph) = — [ [ B{x) 1 / 2 dids{x) + 0{h 1 / 16 ) 

for 6q < b. The proofs of the upper and the lower bound in (|3.12|) are similar 
and we only give the latter one. 

Let N, S be as in Subsection 13.41 and let S n G [s n -i, s n ] be arbitrary with 
B n := B(S n , 0). Note that for all sufficiently small h one has bo — Ch5~ 2 < 
b<b' < B n . Hence, from fl2£3 and Proposition EH for all 5 G (0, 5/2] 



1 N 



n=l 



x fu+ (B-^bo - CM" 2 )) - v_ (b~ x % - Ch5~ 



12 



2 

Now we use (j22J), choose 5 = h 7/ie and recall that N = 0(h~ 3/8 ) to get 

JV 

2 IT 



1 

h}/ 2 N{hb ,P h ) > Y^SBl' 2 (v + (B-%) -^(B-%)) -c^ 1 / 16 . 



n=l 



The main term on the RHS is a Riemannian sum. Recalling that the S n 
were arbitrary, we finally arrive at 

h l / 2 N{hb ,P h ) 

= —! I B{x) 1 / 2 d£ds(x)- ci/i 1 / 16 . 

This is the lower bound of (|3. 12|) . In a similar fashion one can establish the 
upper bound, which concludes the proof of Theorem ll.il 
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4. Proof of Theorem 11.21 

4.1. The model operator on a half-cylinder. For parameters S, T > 0, 
k£| satisfying 

2\k\T < 1, 

we denote by M'^ ,T,K the self-adjoint operator in .^((O, S) x (0,T),a K dsdt) 
associated with the quadratic form 

m S h ' T ' K [u]:= / (a- 2 \(hD s +t- Kt 2 /2)u\ 2 + \hD t u\ 2 ) a K dsdt, 
Jo Jo 

V[fh S h ' T ' K ] := {u G H\(0, S) x (0, T)) : u(.,T) = 0, u(0, .) = u(S, .)}. 

Recall that the function a K was defined in (|2.5|) . (We emphasize once more 
that here and in the next subsection, k will be a constant and not the 
curvature.) 

Moreover, recall the constant C\ from (|1.5|) . The goal of this subsection 
is to prove 

Proposition 4.1. Let D > 0. T/ien i/iere exist C, e > suc/i t/iai /or a// 
N < D, \k \ < D, S > 0, < h < e and e^Vh] logh\ < T < e/i 1 / 4 one has 



< a 



For the analysis of the operators M^' T ' K we begin as in the proof of 
Lemma lM.21 By separation of variables and a dilation r = /i _1 / 2 i we obtain 
the unitary equivalence 

(4.1) Mf T ' K /i^©Af(2vrn/i 1/2 5- 1 ,/i 1/2 K,/i- 1 / 2 T). 

Here we define, for parameters 

(£K, a €[-1,1], L>1, 2|a|L<l, 

the self-adjoint operator M(£,a,L) in the Hilbert space ^((0, L), a K dr) by 
the quadratic form 

m(£, a, L)[/] := T (|/'| 2 + a~ 2 (£ + r - ar 2 /2) 2 |/| 2 ) a a dr, 
•/ o 

2>[m(£,a,L)] := {/ G ^(0,L) : /(L) = 0}. 

The proof of Proposition 14.11 relies on the following two results, which we 
take from HeMo (Section 11), see also [FoHelj (Lemma 5.4). We shall 
denote the eigenvalues of a self-adjoint and lower semibounded operator T 
with compact resolvent by (J-i(T) < (J-2(T) < taking multiplicities into 
account. 
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Lemma 4.2. Let D > 0. Then there exists a C > such that if \a\L 2 < D 
then for all j E N 

M M (C, a, L)) - ^(M(C, 0, L))\ < C\a\L 2 (l + ^(M(£, 0, L))). 

The previous lemma follows easily by comparing the corresponding qua- 
dratic forms. Using an explicit trial function one can show 

Lemma 4.3. Let D > 0. Then there exists a C > such that for all 
\£, — (,o\ < D , L > C there exists a A G <r(M(£, a, L)) w't/i 

|A-e -3Ci|£ |(£-&) 2 + Cia| 

< C (|e - Co| 3 + |a||C - Co| + a 2 + e"^ 
We will use the following consequence of the two preceding lemmas. 

Corollary 4.4. Let (3 > 0. T/ien i/iere exist e,8,C > swe/t £/iat /or a// 
C G M, |a|L 2 < e, L > er 1 one has: 

(1) // |£-6I ^C" 1 fften 

| m (M(e, a, L)) - 6 - 3Ci|£o|(£ - Co) 2 + C x a\ 

< C ( |e - Col 3 + |a||C - Col + a 2 + e- L / c ' 



(4.2) 



and 

/i 2 (M(£,a,L)) >Q + 5. 
(2) K - Col > «ien 

m (M(£,a,L)) >e + 5. 

The parameter /9 is introduced for technical reasons which will become 
clear in the proof of Proposition 14.11 

Proof of Corollary \4-4\ By Lemma 14.31 there exist e, C > such that if |C — 
Col < /?C _1 , |a|L 2 < e and L > e~ l then M(C, a, L) has an eigenvalue 
A below ^(6o + ®i) satisfying (|4.2|) with the constant C. (We could of 
course choose e = C _1 or e = /3C _1 , but later it will be useful to keep them 
separated.) Note that by the variational principle we have 

(4.3) H (M(£,0,L)) > H {M{i)), 

where M(£) denotes the operator studied in Subsection 12.11 (corresponding 
to a = and L = oo). Hence by Lemma 14.21 we find that, after decreasing 
e if necessary, one has for all |£ — Col < /?C _1 

M 2 (M(C,a,L)) > (1 -ci|a|L 2 )^ 2 (M(C,0,L)) -ci|a|L 2 

> (1 - ci|a|L 2 )/x 2 (M(C)) - ci|a|L 2 

> (1 -c 1 \a\L 2 )Q 1 — c\\a\L 

> (9 + 6i)/2. 

In particular, it follows that A = (j,i(M(£, a, L)). This finishes the proof of 
the first part of the corollary (with 5 < ^(@\ — Go) arbitrary). 
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By the properties of the function [i recalled in Subsection 12.11 it is clear 
that there exists a constant 6 > such that for all |£ — £o| > /3C _1 one has 

M(£) > ©o + 2<5. 

(We can assume that 5 < ^(0o + ©i).) Applying again IJ4.M|) and Lemma l4~2l 
and decreasing e if necessary we find for all |£ — £o| — /3C _1 

Hi(M(£,a,L)) > (1 - ci|a|L 2 )^i(M(^,0,L)) - ci|a|L 2 

> (1 - ci|a|L 2 )/u(^) - ci|a|L 2 

> (1 - ci|a|L 2 )(9 + 25) - ci|a|L 2 

> e + s. 

This finishes the proof of the second part of the corollary. □ 

Everything is now in place for the 

Proof of Proposition ^. 1\ We keep the notation e, <5, C for the constants from 
Corollary 14.41 corresponding to 

(4.4) := 3Ci|&|/2, 

and we will assume that 

CVh\ log h\ < T < ^/TjD h 1/A . 

Then there exists a ho > such that for all < h < ho and all |«o| < 
D one has h~ l / 2 \K\T 2 < e, hr x l 2 T > e" 1 and h 1/2 CiK < 5. Therefore 
Corollary 14.41 and the decomposition (|4.1|) imply that 

iv(/ie + /i 3/2 CiKo,M^ TA ) 

= J{neZ: fH{M(27mh^ 2 S- 1 , h 1 / 2 ^ h~ l / 2 T)) < 6 + /i 1/2 Cik }. 

Noting that hn 2 + e ~ T/c ^ < h(D 2 + 1) we find from Corollary E3 the 
estimates 

|{iiEZ: p + (|2 7 rn5- 1 - /i~ 1/2 £o|) < 0} 
< iV(/i6 + /i 3/2 Ci^o, Mh' T ' K ) 
"" <|{n£Z: ^vrnS" 1 - /i~ 1/2 £ | < PCT 1 ^ 1 ' 2 , 

p-C^miS -1 - /i- 1/2 £ |) < 0} 

where 

P±,fc(y) := 3d|£o|y 2 - Ci/T 1 ^* + «q) ± C(/ l 1 / 2 y 3 + + D 2 + 1). 

The assertion will follow from the properties of these polynomials which we 
will discuss now briefly. 

We start with p +>h . If k + k < ^(Z) 2 + l)/! 1 / 2 one has p +)h (y) > 
for all y > and we define := 0. On the other hand, if k + kq > 
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C^ 1 C(D 2 + l)h 1 / 2 , one checks that there is a unique zero y +: h G (0, oo), and 
that this satisfies 

(4.6) y +A = (3|e |)" 1/2 (^ + Kaffir 1 /* + 0{l) 

as h — > 0, where 0(1) is uniform in k, ko (varying in a bounded set). Hence 
in any case, we obtain from (|4.5|) the lower bound 

N{h& + h 3 / 2 C lKo , M^ ,T,K ) > JJ{n G Z : ^nS" 1 - /T 1/2 £ | < y+,4 

= /»-V4_|^ (f6 + ISo) i/2 + 

Now we turn to the polynomial p_ First we note that there is a 
zero y~ t h ~ 3Ci|£o|C -1 ^ -1 ^ 2 ) which however does not lie in the interval 
[O^C^h- 1 / 2 ) if h is sufficiently small. (This is the reason for the choice of 
P in (gai).) If k + k < -C^ l C(D 2 + l)/i 1/2 , y_ ih is the only zero in (0, oo) 
and we set y^^ := 0. On the other hand, if k + ko > —C± l C{D 2 + l)/! 1 / 2 , 
one checks that there is a unique zero y~h G [0,y /J, and that this has the 
same expansion as in (14. 6j) . In both cases, we obtain the upper bound 

iV(/»e + h 3/2 dK ,M^ T ' K ) < (t{n G Z : ^nS" 1 - /T 1/2 £ | < !/-,/»} 

= ^ 1/4 ^^(^ + ^o)f + 0(l). 

This proves the assertion. □ 

4.2. The model operator on a Dirichlet strip. We fix S, T, k as in the 

previous subsection and consider the operator M h ' ' obtained from M h ' ' 
by imposing additional Dirichlet boundary conditions at s G {0,5}. More 
precisely, M^' T ' K is the self-adjoint operator in ^((0, S) x (0,T),a K dsdt) 
associated with the quadratic form 

m ^ K [ u ]:= / ( a ~ 2 \(hD s + t- Kt 2 /2)u\ 2 + \hD t u\ 2 )a K dsdt, 
Jo Jo 

V[m S h ' T ' K ] := {u G ^((0, S) x (0, T)) : u(.,T) = u(0, .) = u(S, .) = 0}. 
With an argument similar to that in Subsection 13.21 one proves 

Proposition 4.5. There exists a C > such that for all S,T > 0, k£I 
with 2|k|T < 1 and all A 6l and 5 G (0, S/2] one has 

-N(X - C/i 2 <T 2 , A? 2(5 ~ 5) ' T ' K ) < N(\,M%' T ' K ) < N(X,M^' T ' K ). 

4.3. Estimates near the boundary. Similarly as in Subsection l3.3l we will 

now approximate the quadratic form locally near the boundary but, the 
magnetic field now being constant, this can be done with a higher precision. 
In particular we will see the curvature of the boundary appear. Again we 
shall use the notation from Subsection 12.21 and, in contrast to Subsection 
13. 3( it will be important now to keep k arbitrary. 
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We will assume in this subsection that B = 1. Then we can choose the 
magnetic vector potential in the following way. 

Lemma 4.6. There exists a constant C > such that for all S G (0,L), 
S G [0,5] there exists a <j)<E C 2 ([0, S] x [0,i ]) such that 

A(s, t) - V(p{s, t) = (-t + kt 2 /2 + (3(s, t),0) T , (s, t) G [0, S] x [0, to], 

where k := k(S) and for any < T < to 

(4.7) sup \(3(s,t)\ < CST 2 . 

(s,t)e[a,s]x[o,T] 

Indeed, one can take (p(s, t) := J Q * d s A 2 (s, t') dt'—J Q s A\(s' , 0) ds' and recall 
that d s A 2 (s,t) - <9 t ii(s,t) = 1 - t«(s). 

Before stating the next result we recall the definition of v from (|2.6|) and 
of the quadratic form mf i ,T ' K from Subsection 14.21 

Lemma 4.7. Let D > 0. Then there exists a constant C > such that for 
all S G (0,L), S G [0, S], T G (0,t ] > ^v 7 ^ and /or a// u G 2% ft ] 

smc/i i/iai the corresponding v satisfies Q3.6JI one /ias 

Here k := k(S') and <j) is the function from Lemma \4-b\ Moreover, \\ ■ \\a 
denotes the norm in L2((0,S) x (0,T), a^dsdi). 

Proof. The proof is rather similar to that of Lemma 13.51 so we will only 
sketch the major steps. Writing w := e~ l( ^l h v and taking (|2.7|) into account 
we decompose 



q h [u] - m^ T ' K [w] = Lx + L 2 + L 3 



where we define 

"t 



h-= J (cT 2 \{hD s - ii)w| 2 + \{hD t - A 2 )v\ 2 ^ a k dsdt 



r-C rto 

h 2 / W k \v\ 2 a~ K dsdt, 



o Jo 
h 2 f c 



k - k)\v(.,0)\ 2 ds. 



To treat the terms L 2 and J3 we use that \n — k\ < c\S on the support of w. 
This leads to the estimates 

\L 2 \<c 2 h 2 (S + T)\\w\\l 

\h\ < e\\hD t w\\% + c z e- l h 2 S 2 \\w\\l. 
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for any e > 0. To take care of I\ we note that by Lemma 14.61 one has 

h= f r\a-\hD s + t-Rt 2 /2-l3)w\ 2 
Jo Jo 

-a^ 2 \(hD s + t - kt 2 /2)w\ 2 ) a k dsdt. 

Using the estimate on (3 from 1)4.7)1 and that \a~ 2 a\ — 1| < c^ST on the 
support of w we easily find, for all e > 0, 

rC /-to 

\h\ < (e + c 5 ST) / a^ 2 \(hD s + t- kt 2 /2)w\ 2 a k dsdt 
Jo Jo 

+ C5 (i + e- 1 )s 2 r 4 H|i. 

The assertion then follows by choosing e = ST and recalling T > D\f~h. □ 

4.4. Bracketing. Now we estimate N(hQo + h?l 2 C\Ko, Ph) by the spectral 
counting functions of the operators with constant curvature from Subsec- 
tion ^3 Again we assume that B = 1. 

For N G N we define S and s n as in 1)3.1 In contrast to Subsection 13.41 
we will not yet specify the value of N but postpone this to the next subsec- 
tion. 

Proposition 4.8. Let kq G R. Under the above hypotheses there exists 
a constant C > such that for any < h < C" 1 , 5 G (0, S/2], T with 
Ch^ 2 \logh\ < T < C^h 1 / 4 and S n G [s n _i,s n ], n = l,...,JV, one /ias 



- £ N(h@o + /i 3/2 Ci«o - C(ST 3 + ^ 2 <T 2 ), M ft 2(S+J) ' T '^) 

2 n=l 

<iv(/ieo + /i 3/2 Ci K o,n) 

AT 

< N(hO + /i 3/2 Cik + C(ST 3 + h 2 5~ 2 ), M^ +2S ' T ' Kn ) 

n=l 

where k n := n(S n ). 

The proof is similar to that of Proposition 13.61 where however Proposi- 
tion l3.3l and Lemma f3.4l are to be replaced by Proposition ^ . 51 and Lemma f4.6l 

respectively. We omit the details. 

4.5. Proof of Theorem 11.21 Replacing h by h/B we may assume that 
B = 1. We will show that 

i r c 

h l ' A N{hQ + h^CtKo, P h ) = — == / (k(s) + kq) 1 / 2 ds 
(4.8) vrVS Jo 

+ o e {h 1 / w - e ) 

for any e > 0. As in Subsection 13.51 we give the proof of the lower bound 
only. 
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Let N, S be as in Subsection 14.41 and let S n S [s n _i,s n ] be arbitrary 
with k n := n(S n ). From Proposition 14.11 and Proposition 14.81 we get for all 
5 G (0,5/2], Ch l l 2 \ \ogh\ <T< C^h 1 ^ (with C as in Proposition Q|) 

(4.9) 

h^NihQo + h^d^Ph) 

N 

> (Trv^j&l)- 1 ^ - s) + k - Cl (h- 3 / 2 ST 3 + hV 2 r 2 ))¥ 2 

n=l 

- c 2 Nh 1 /\ 
Now we use the estimate 

(k n + K -c 1 (h- 3 / 2 ST 3 + h 1 / 2 6- 2 )) 1 / 2 

> {kn + Ko) y 2 _ c }/2 ( ^-3/4 5 l/2 T 3/2 + fcl/^-l^ 

If we assume that 8 > h 1 / 4 we easily deduce from (|4.9j) that 
(4.10) 

N 
n=l 

- c^h-^S 1 ' 2 ^' 2 + h l 'H- 1 + 5S- 1 ). 

Now choose T = h l l 2 ~ p with 0<p<l/4. A calculation shows that the 
second term on the RHS of (|4.1()|) is minimal for the choice 

£ = /jl/8+3/V2 g = h 3/16+3p/4 

and given by C4/i 1//16_3P//4 . The first term on the RHS of 1)4.10(1 is a Rie- 
mannian sum. Recalling that the S n were arbitrary, we finally arrive at the 
lower bound of 1)4.8)1 . The upper bound can be established similarly, which 
concludes the proof of Theorem ll.2l 
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